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Abstract 

We study derealization transition in a one- dimensional electron system with quenched 
disorder by using supersymmetric (SUSY) methods. Especially we focus on effects of 
nonlocal correlation of disorder, for most of studies given so far considered ^-function 
type white noise disorder. We obtain wave function of the "lowest-energy" state which 
dominates partition function in the limit of large system size. Density of states is cal- 
culated in the scaling region. The result shows that derealization transition is stable 
against nonlocal short-ranged correlation of disorder. Especially states near the band 
center are enhanced by the correlation of disorder which partially suppresses random 
fluctuation of disorder. Physical picture of the localization and the derealization 
transition is discussed. 
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1 Introduction 



Quenched disorder plays important roles in various physical phenomena. Anderson 
localization is one of these examples Especially in two and lower dimensions almost 
all states are localized, and extended delocalized states are isolated points in physical 
parameter region. The transition between quantum Hall plateaus is such an example. 

For the transition between integer quantum Hall plateaus, useful model, named 
network model, has been proposed^], which incorporates effects of localization and 
quantum tunneling in a strong magnetic field. Numerical studies indicate that the 
network model belongs to the same universality class of the transition between integer 
quantum Hall plateaus. However it is rather difficult to solve the network model 
analytically because there exists no controllable parameter for perturbative expansion 
and also because of nature of quenched disorder itself. Compared with the two- 
dimensional (2D) systems, one dimensional (ID) systems are more tractable. These 
ID systems include Dyson's study on random strength harmonic springs 0, random 
Ising model by McCoy and Wuj|], and more recently random exchange spin chains Q 
and random hopping tight-binding (RHTB) model|6|, [7|, ||]. 

Recently supersymmetric (SUSY) methods appear useful for handling the quenched 
disorder. SUSY methods are applied to the network model||, the ID RHTB model0 
etc. In this paper we shall revisit the ID RHTB model by applying the SUSY meth- 
ods. Especially we shall consider nonlocally correlated quenched disorder and study 
stability of the derealization transition which exists in the case of the 5-function-type 
white noise disorder. The model contains two parameters which control magnitude 
of fluctuation and correlation length of disorder. We expect that we can get detailed 
physical picture of (de) localization transition from calculations of density of states, 
Green's functions, etc. as a function of the above parameters. Another motivation of 
the present work is rather technical, i.e., we show how to use the SUSY methods for 
nonlocally correlated disorder systems. 
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Hamiltonian of the RHTB model is given by 

HRHTB = -XX( C n C n+l + C Ul C n), (1-1) 
n 

where c n is annihilation operator of spinless fermion at site n and t n 's are random 
hopping parameters. As we shall see, staggered part of fluctuation of £ n 's plays an 
important role at low energies and generates energy gap for most of fermion modes. 
Extended excitations are located near the band center, and they are described by a 
Dirac fermion with randomly varying mass. 

In the SUSY methods, bosonic variables are introduced. In terms of them, density 
of states, Green's functions, etc are expressed in a compact way. Average over random 
variables of nonlocal correlation can be performed by introducing another bosonic 
variable, and because of that there appear couplings between the fermion and its 
bosonic SUSY partner. The model is reduced to a quantum mechanical system of 
one fermionic and two bosonic variables. The lowest-energy state which dominates 
the partition function is annihilated by SUSY charges. We obtain wave function of 
the lowest-energy state and calculate density of states near the band center. 

This paper is organized as follows; In Sect. 2, we shall review the RHTB model, 
its continuum limit and SUSY methods which are applied for nonlocally correlated 
disorder. Most of our notations follow those by Balents and Fisher in Ref.J?]]. In 
Sect. 3, transfer Hamiltonian is obtained regarding the spatial coordinate as time. 
Then the system is reduced to a quantum mechanical system. In Sect. 4, we shall 
obtain an equation of motion of the ground state, which dominates the partition 
function and density of states. Section 5 is devoted for solution to equation of motion, 
and in Sect. 6 by using the ground state solution we shall calculate density of states 
of fermions. Section 7 is devoted for discussion. 

In this paper, we calculate density of states of fermions. We shall report Green's 
functions in a future publication [|K|. These calculations are quite useful for under- 
standing localization and the derealization transition. 
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2 Continuum limit and SUSY 



The random hopping parameters fluctuate as 

t n = t + 5t n , (2.1) 

where t is some finite constant and we normalize as t = ^- and a\ is the lattice spacing. 
We assume that r.m.s. of random variables 6t n is sufficiently small compared with the 
band width t. Continuum limit of the system ( |Llp is then easily obtained. We shall 
focus on modes near the band center because derealization transition occurs there. 
Excitations near the band center of fermions are described by smoothly varying right 
and left-moving fermion fields ipR and ipL, 

c n = j kFH Mx) + e~ lkFn Mx), k F = |, (2.2) 

where continuous coordinate x = na\. 

The most relevant part of the random hopping parameters St n is their staggered 
part, 5t n ~ (— l) n m(x). Then in terms of ipn(x) and we obtain continuum 

limit of the Hamiltonian (|1 . 1| ) as 

i>Rid x ipR ~ tpL id xipL - im(x)(i) R 4) L - 4>{iPr)\ , 

h = -ta z d x + m(x)a y , V = (Ml)'- (2.3) 
The random mass m(x) is decomposed into a uniform and random piece as 

m(x) = m + <f>(x), (2.4) 

where [<fi] ens = and 

[(j){x)(j){y)} ens = — exp(-|x-y|/A), (2.5) 
with positive parameters g and A. It is easily verified 



H c = - /dx 



ydx [<f)(x)<f)(y)] ens = g, 
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and 

[4>(x)(f>(y)} ens -»• g8(x - y) as A -> 0. 

It is obvious that the parameter g controls magnitude of fluctuation and A is the 
correlation length of disorder. 

Single fermion Green's function at energy u is defined as 

G a p{x,y;iu) = (x,a\ \y,P), a,P = R,L, (2.6) 

where \x, a) is a normalized position eigenstate of fermion at x and chirality a. By 
functional integral, 

G aP (x,y;iu) = %{^ a {x)^p{y))^ 



S^p = /dx ip(ih + uj)ip, 

Z$ = [ Vij)Vij)e- s + . (2.7) 



Ensemble averaged Green's function is obtained from ( |2.6j ) as 



G af3 (x-y;iuj)=\(x,a\ 1 . \y,(3)} , (2.8) 

L h — IUO J ens 

where the ensemble average is taken with respect to <f)(x) in ( |2.3| ) and ( |2.4| ) according 
to ( |2.5D . To take the ensemble average seems formidable. However by introducing 
SUSY partner, Green's function G a p(x — y; iuS) etc can be expressed by a functional 
integral in a compact way. Point is that integration over the bosonic SUSY part- 
ner cancels the fermionic determinant and normalization of the fermionic states is 
automatically taken into account. Introducing the bosonic superpartner £, 

G af3 (x - y; iu) = i(^ a (x)^)p(y))s, (2.9) 

where 



(2.10) 



with 

S = Jdx [$(ih + + £(ih + (2.11) 

The ensemble average in ( 2.10| ) can be performed in the following way. We first 
notice identity such as 

(-X 2 d x 2 + 1)^ exp(-|x - y\/X) = 8{x - y). (2.12) 

Then the ensemble average can be converted into the functional integral form, 

[0(aO0(3/)...U = j V<P(x')(<P(x)<P(y)..^M-SM*')])i (2-13) 
S*[<f>{x)) = / dx ^(p(x)(-^d x 2 + l)4>(x). (2.14) 
From ( 2.10Q , (|2.11|) and (|2.13 ), the expectation value of operator A is given by 

(A) s = J Vi)Vi)ViViV<\> A exp(-(3 + S^)). (2.15) 

The above total action is obviously SUSY under ip <-> £ and has such a form that 
the SUSY partners ip and £ couple to the "dynamical" real scalar field (f) which is 
"singlet" under SUSY transformation. 

In the following section, we shall obtain a transfer Hamiltonian by regarding the 
spatial coordinate x as time. The system reduces to a quantum mechanical system 
with the two bosonic and one fermionic variables. 



3 Transfer Hamiltonian 

In this section, we shall obtain a transfer Hamiltonian by regarding the spatial co- 
ordinate x as time in the functional integral representation Q2.15 ). Then the system 



reduces to a quantum mechanical system. In Ref. [/J, the following canonical creation 
and annihilation operators are introduced corresponding to the functional integral 
variables, 

i> R ^F h fa^-Fi, ?Pl->fI 

&-fl T , U^B\, h^~B\. (3.1) 
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It is useful to define fermionic and bosonic spin operators, 



j b = \bbb, 



where 



B = B ] 



a" 



(3.2) 



(3.3) 



It is proved that Jf and J b satisfy SU{2) and SU(1, 1) algebras, respectively, and 
they commute with SUSY charges Q and Q, 



Q = BF, Q = F^B. 



(3.4) 



Transfer Hamiltonian of F a and B a (a =t, I) part is written in terms of the spin 
operators. 

Transfer Hamiltonian of <fi is also obtained from ( |2.14|) . The system of (f) is nothing 
but a simple harmonic oscillator linearly coupled with the SUSY spin J = Jf + J b . 
In terms of the spin operators and the canonical boson operators of the harmonic 
oscillator a, a) which correspond to 0, Hamiltonian of the system is given as, 



H = 2uJ z + 2m J x + \%J x {a + a)) + \{<Ja + \). 

V A A 2 



(3.5) 



Fermionic states of F a are specified by representations of SU(2). Similarly bosonic 
states of B a form multiplet of irreducible representations of SU(1, 1), which are spec- 
ified by total spin 



J 2 = {N 2 B + 2N B )/4, N B = BB = B}B^ - B\B b 



(3.6) 



and ^-component of spin J z , i.e. 

J 2 \jn) 
r\jn) 



l + |2j + l| 



+ n 



(3.7) 
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where it should be remarked here that the total spin takes half integers, namely, 
j=0,±l/2,±l,--. 

In Ref. |7j , structure of the above quantum mechanical states of the SUSY partners 
is studied in detail and it is shown that SUSY-invarinat states are explicitly given as 
follows; 

f i \\-l/2,n)®\l) F +\-l/2,n-l)®\ n>0 
n) = { V (3.8) 
[ | - 1/2,0) <g) | i) F , n = 

where \u)f is the fermionic state of spin a. 



4 The ground state 



From the functional-integral representation (|2.15| ), it is obvious that the partition 
function Z = Ti((—l) N fe~ LH ), where Nf = fermion number and L = system size, 
acquires contribution only from the scalar field <fi with the action because fermionic 
and bosonic determinants cancel with each other for the SUSY partners. Then in the 
large system size limit L — > oo, Z = e~ L ^ 2X . The "ground state" which satisfies 

#|0> = ^|0>, (4.1) 

or 

(2ujJ z + 2m J x + \j\j x (a + a f ) + -^a) |0) = 0, (4.2) 

dominates the partition function and the density of states of fermions which we 
shall calculate in this paper. In Eq. (|4.2| ), we have rescaled the parameters as uj — > 



2gco>, tuq — > 2gmo and A — > ^. This rescaling leads the general Hamiltonian to the 
one of g — |. After calculation, we shall scale back the above parameters and recover 
dependence on g. 

The SUSY vacuum should be invariant under SUSY transformation and it must 
be annihilated by both Q and Q, and therefore it is given by a sum of a direct product 
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of \n)o in Eq.( |3.8| ) and \rn)H, which is the eigenstates of the harmonic oscillator, i.e. ; 

1°) = H <f>n,m\n) \m) H , (4.3) 

n,m 

where = ^/=r ( a ^) m |0)-ff with the "vacuum" of the harmonic oscillator |0)#. 

Norm of the state ( [4.3|) is given by[| 

(0|0)=El^o, m | 2 . (4.4) 

m 

Let us insert Q into O, 

Y,4>n,mhujn\n)o\m) H + m ((n + l)|n + l) - (n - l)|n - l) )|m)^ 



+ y ^((n + l)|n + 1) - (n - l)\n - l) )(v / m + l|m + 1)^ + v^|m - 

+ ^m|n) |m)//^ = 0. (4.5) 

Rearranging the indecies n and m of the subscripts of 4>n,m an d the labels of states, 
we obtain equation of <p n ,m for the vacuum, 

(2un + ^m)0„ im + m • n ■ {<f) n -i,m ~ 4>n+i,m) 



+ y— n(VmTT(i-l, m+ l - <fi n+ l,m + l) + V / ^l(0n-l, m -l - 0n+l,m-l)) = 0. 

(4.6) 

Equation ( ^4.6| ) for n = can be solved easily and the solution is 

0, m^O, 

0o, m = (4.7) 

1, m = 0, 

where the normalization comes from the condition (0 1 0) = 1. 

By solving Eq. Q4.6j ), we can calculate the density of states of fermions. However, 

it is not easy to obtain solution for arbitrary A. Then in the following section, we 

shall solve Eq. (|4.6|) by perturbative calculation in powers of A, for A = corresponds 

to the case of 5-function type white noise. 

1 As shown in Ref(7j], right and left states must be distinguished in the present system because of 
the non-hermiticy of the Hamiltonian. Then only (j) n =o,m contributes to the norm. 
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5 Solution 



In this section we shall obtain wave function of the vacuum state by solving (|4.2| ). We 
consider the case in which A is small. The limit A — > corresponds to the short-range 
limit which is discussed in Ref.[[F|. 

In order to see the correspondence between the transfer Hamiltonian of the present 
system ( |3.5| ) and that of the short-range limit studied in Ref. [|7| , 

H = 2uJ z + 2m Q J x - 2(J X ) 2 , (5.1) 

we rewrite the Hamiltonian ( |3.5|) as follows, 

H = H + \(a^ + V2XJ x )(a + v / 2\J x ) + ^-. (5.2) 
A 2 A 



From (|5.1| ), it is easily seen that the second term in Eq. ( |5.2j ) should vanish in the 
A^O. 

The ground state |0) of the transfer Hamiltonian (|3.5|) will be obtain by the 
perturbative calculation from the state lO)*- -*, which is annihilated by the operator 
(a + V2XJ x ), 

(a + V2\J X )\0)W = 0. (5.3) 

Decomposing the state |0)^ into the states \n)o\m)H, namely, \0)^ = 4>n) m \ n )o\' m )H, 
Eq.( |5.3| ) leads to equations of </>^, 

v^TT^L+i + V - <j> { nl l>m ) = 0. (5.4) 



Equation (|5.4|) shows that, in the limit A — > 0, 4>^ m = for m > 1 and therefore 



4>nj) = 4>ni where <p n is the wave function for A — >• defined in Ref.0, i.e., 

ffoE^n|n)o = 0. (5.5) 

n 

We shall consider higher-order corrections of finite A. Equation (|5.4|) can be easily 
solved as 

1 (2A)f(nA n )">!S, (5.6) 
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where A n is the difference operator with respects to the index n, i.e., A n (f)£ 

-*C0) _^(o) 

-, and this tells us the order of 



'n+l.m 'ti — l,m 



2 , «^j.j.vj. uij.x>j uij.^ vyj. ^rn^mi 

At. (5.7) 



Now let us turn to the estimation of the difference between the states |0) and \0)^°\ 
We decompose the former as 

\Q) = ^ { nl + ^n,m)\n)o\m) Hl (5.8) 



where 5(fi n , m comes from difference between |0) and |0)(°). Substituting (|5.8|) into ( flT] 
we obtain relation between 4>^ m an d ^0n,mi 



^(tfoA + # A <ty»,nO|n)o|m)fl = 0, (5.9) 

n,m 

where H X = H - i Noting # A |0) (0) = #o|0) (0) , we rewrite flp|) as follows, 
E (Wgi. + (#o + 2(JT) 2 + y)50n, m 

n,m 

+]l^J x (VmS(pn,rn-i + V m + 150„, m+ i)) |n) |m}# = 0. (5.10) 
For each m, the terms should be the same order of A, i.e., 



H 0<Pnm ~ jHn,m ~ ~ ^=(^0n,m-l + ^n,m+l)- (5.11) 



From the fact that ifo does not depend on A and ( |5.7|) , we obtain 

Hn, m ~ A(t+i), (5.12) 



which is of higher order of A than 4>^ m f° r the same m. From the above result, we 
can estimate order of the term depending on A in ( |5.2j ). From Eqs.( |5.7|) , (|5.12 ) and 
the following equation, 



a T + V2XJ x )(a + V2XJ X )\0) 



E( m ^n,m + v / 2AnA n ( v /m50 njm _i + Vm + i5<f} n , m +i) 

n,m 

+2\(nA n ) 2 6(f) ntm )\n) \m) H , (5.13) 
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we find that the above term is O(X^). Then it is clear that our estimation is consistent 
with the expectation that in the limit A — > the present system becomes that studied 
in RefJ]. 

In order to calculate the density of states, we have to solve Eq.(O). Decomposing 
it for each powers of A might be the easiest way to find 4> n ,m- 4>no contains higher-order 
term of A than n . For m = 0, Eq. (|5.10|) reads 

£ (#o(<® - 0n) + (Ho + 2(jP) 2 )^n,o + ^J'V^+IS^ |n)o = 0, (5.14) 

where we use the fact that Ho annihilates the vacuum for A = 0, and one can find 
that <f$ Q - (f) n is 0(A). Then from we can see that 0^ contains a term of 

0(Air +1 ) as the next-leading-order term, and therefore we can decompose it as 0^ = 
0n,m + A0„ jm , and 4> n ,m and 4> n ,m are both C(A^). We notice the order of A in Eq.( |5.6|) 
and find 

n , m = -L=(2A)T(nA n ) m n . (5.15) 



We insert the above ^nln Eq.( |5.10| ) and extract the terms in the order of 0(X 2 +1 ) 
and C(Af ) ; 

(^H 4> n!m + (H + 2(J x ) 2 )b<\)^ m + ^J x Vm + l<ty n , m+ i j |n) = 0, 

^Ho4>n,m + y^0n,m + \^J* V^<t>n,m-l \ \n)o = 0- (5.16) 

By inserting the second equation of (|5.16|) into the first one, we obtain 



/ 2A 

Ho(\(fin,m + 5(j)n,m)\n)o = J — — ^J x Ho4> n ,m+\\n)o 

1-) V I -L ,n 



2A 

' T E^^O^0n, m |n) O , (5.17) 



m + 

where we have used Eq. (|5.4|) in 0(\~~5~), namely, 

Vra + 1^0 n>m+ i|n) o + \p2\J x Y<p n , m \n)o = 0. 
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We shall solve Eq.( |5.17|) from now on. 



Let us consider the m = case first, where n> o equals n . Appendix C in Ref.fTj 
gives us the solution (fi n . We review it here for completeness. We recall that J2 n 4>n\n)o 
is the vacuum state of the transfer Hamiltonian with A = 0, and thus <p n satisfies the 
following equation, 

hum - 2m nA n - 2(nA n ) 2 ) n = 0. (5.18) 
If <p(n, mo) is a solution to Eq. (|5.18 ), then so is (— l) n <ft(n, —mo), and thus general 



solution can be written as 

<Pn = Ci0(n, m ) + c 2 (-l) n 0(n, -m ), (5.19) 

with constants c% and ci- Since <p n = 1 and <p n = (—1)" are solutions of the above 
equation for u = uiq = 0, <p(n, mo) for small to and mo is a slowly varying function of 
n. This suggests that <p(n, mo) is a "continuum" function of n and that enables us to 
replace discrete differences with derivatives. The differential equation for (f>(n, mo) is 
thus obtained as, 

(2un - 2m nd n - 2(nd n f) <p{n, m ) = 0. (5.20) 

In order to solve the above equation, we embed the integer n into complex number 
and use the Laplace transformation, 



(j)(n,m ) = / dt e~ nt (j){t). (5.21) 
J o 

By inserting the above Eq.( |5.21|) into Eq. (|5.20 ), we have 



u<j>(t) - (1 - m Q )t(j>(t) - t 2 d t (f)(t) = 0, (5.22) 

with the boundary condition <f)(t) —>■ in the limit of t — > or t — > oo. This boundary 
condition comes from the requirement that <p(n,mo) tends to vanish for large n. We 
find the solution 

roo 

4>(n,m ) =a dt ^o-i e -nt- T) ( 5 _ 23 ^ 
Jo 
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where ao is a constant, and the above solution can be written in terms of the modified 
Bessel function 

K„(z) = K_ u (z) = -Qy / dt e-*-**-*- 1 (5.24) 



2 y 2' Jo 

as 



4>(n,m ) = 2a (-) 2 K mo (2Vum). (5.25) 



m' 

In order to obtain the constants CiOq and C2CI0, we use the boundary condition |7|, 



0o ~ ci0(l, m ) + c 2 0(l, -m ) ~ 1, (5.26) 
A n 0i ~ ci<9 n 0(l, m ) - c 2 d n (f)(l, -mo) ~ 0. (5.27) 



Please note that the modified Bessel functions have a singularity at n = 0. Eq. ( [5.18 



is, therefore, ill-defined for n = 0, but we can expect the "genuine" 0(0, mo) in O 
has almost the same value of 0(1, mo), since the function 0(n, m ) is assumed to be 
a slowly varying function. For small n, we can approximate 0(n, mo) as 

0(n,m o ) 

/■DO 

= a / dt f"o-i e -nt- T 
Jo 

- a ° (/ dtr°-V^+ / dtt^-V*) 
-(r(m ) + (^) mo r(-m )) 



n m ° Vo 
a 



n m 
Qq ( 
m 



(5.28) 

which coincides with the result obtained in the hard- wall approximation ]7]. Using 
this, Eqs.flHD and (113) S ive 



^(1 - ^) + - = 1, (5.29) 

Cia = C2ao- (5.30) 

Thus the constants are 

moUJ mo 

Cl a = c 2 a = 1 _ uj2mo - (5.31) 
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Now let us return to Eg. ( |5.17| ) . Renumbering the indecies of each terms, we have 

(2un - 2m nA„ - 2(nA n ) 2 )(A0 n , m + 5<f> n>m ) = 
2A 

nA n (2ujn — 2m nA n — 2(nA n ) 2 )nA n nm . (5.32) 



m + 1 

Since 0„ im can be decomposed into the form 

4> n ,m = ci0(n, m, m ) + c 2 (-l) n+m 0(n, m, -m ), (5.33) 

we find that \<fi n ,m + 5(j> n ,m can be also written as 

A0n, m + 50 n , m = A(ciV>(n, m, m ) + C2(-l) n+m ^(n, m, -m )). (5.34) 

As it is for n , we can assume that ip(n,m,mo) is a slowly varying function with 
respect to n. This again allows us to replace the difference operator A n with the 
differential operator d n . Eq.( |5.32|) can be rewritten in the simpler form; 

2 

(2um — 2m nd n — 2(nd n ) 2 )(ip(n, m, m ) H (nd n ) 2 (j)(n,m,mo)) 

m + 1 

2 

2un 2 d n (p{n 1 m 1 mo). (5.35) 



m + 1 

We explicitly solve the above equation for the case m = 0. For n ^ 0, 

(2cu - 2(1 + m Q )<9 n - 2n(d n ) 2 )(ij(n } m ) + 2(nd n ) 2 0(ri, m )) = 

= — 4und n (j)(n, m ), (5.36) 

where ^(n, m ) = ^(n, 0,mo). We again employ Laplace transformation to solve the 
above equation as before, 



i/j(n,m ) + 2(n<9 n ) 2 0(n,m o ) = / dt e -n V(t). (5.37) 

JO 

From ( |5.37| ) and ( |5.21| ), Eq. ( |5.36| ) is rewritten as follows; 

( w _ (i _ mo ) t - t 2 dt)$(t) = 2ud t {t4>{t)) 

= 2a ud t {t mo e-^) 

= 2a uj(m t mo - 1 + ujt m °- 2 )e-^ . (5.38) 
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In general, we can put ip(t) = C(t)e *. Inserting this into Eq.( |5.38| ), we find that 
C{t) satisfies the following equation, 

t 2 d t C(t) + (1 - m )tC(t) = -2ao^(m t mo ~ 1 + ujt m °- 2 ). (5.39) 

A special solution of this equation is C(t) = 2aoUjmot m °~ 2 + aoUJ 2 t m °~ 3 , and thus we 
find 

i/j(n, wo) = -2(nd n ) 2 0(n, m ) + a uj / dt (2m t m °- 2 + ujt m °~ 3 )e- nt -- 

Jo 

= -Aa (nd n ) 2 ((-)^K mo (2V^)) 
n 

fl l-m Tl 2-m 

+2a cu(2TO (-) 2 iT 1 _ mo (2 v ^7) +u(-) 2 K 2 _ mo {2^)). (5.40) 

UJ UJ 

Then we find that the general solution to Eq.( |5.32j ) for m = is given by 

A0n,O + <^0n,O = fHn ~ 2X(nA n ) fa 

fl l-mp 2-mp 

+2Ao;(cia (2TO (-) 2 lfi_ mQ (2vW) + cj(-) 2 i^2- mo ( 2 v^)) 

?7, l+mfl JT, 2+m 

+c 2 a (-l) n (-2m (-)^i ; s:i+ mo (2V^) +w(-) 2 lf 2+mo (2v^))), 

Co 1 UJ 

(5.41) 



where / is a constant which is to be determined by the boundary condition Eq.([4 
For regular n 's, (nA n ) m n = for n = 0, since this term is proportional to n. This 
gives 0o,m = for to > 0, which can be verified by practical calculation. Using fa = 1, 
the condition Eq. (|4.7| ) requires 

A0O,m + 5fa,m = 0, (5.42) 



for arbitrary to. We insert ( |5.41 ) into ( |5.42j ) for the case m = in order to determine 



the constant /. For small n, the modified Bessel function is approximated as 

2{-) l± P 1 K l±mo {2^) - ^ m °-\ (5.43) 
uj 

2{-)- ± ^ L K 2±mo {2V^) ^ u Tm °- 2 - (5.44) 

UJ 
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Then Eq. (|5.42 ) for m = is 



A0o,o + ^0,0 

= fX + X(c l a u}{2m u mo ~ 1 + uuj mo ~ 2 ) + c 2 a uj{-2m uj- ma - 1 + uu~ m "- 2 )) 



m u m ° , 1 - u 2m ° 1 + 



2mo 



Thus 



where 



fX - A - u 2 (2m 

1 4- w 2mo 

/A - 2A™ 2 , + \m Q --—^— = . (5.45) 



4>n,0 = <Pn + A^„, (5.46) 



1 _|_ ^mo 

^ n = (2m 2 , ~ m o-, ^r)0n - 2(nA n ) 2 0n 



1 _ u 2m 



m uj mo 

+2lu-_ 



/ Ti l-m fl 2~m 

2m (-) 2 K x - mo {2yfnuj) + uj(-) 2 K 2 _ m() {2^Jmo) 
+(-l)-(- 2mo (!i)^^ 1+mo (2V^+a;(-)^ir 2+mo (2V^) N ). 



(5.47) 

The other components n)m for m > 1 can be calculated in a similar way to the above. 
Calculations are summarized in Appendix A. 

6 The density of states 

Now we can calculate the density of states. The density of states is obtained from 

G af3 as 

pie) = - limlm£(0;e + ^), (6.1) 
where the averaged Green's function Q (x; iuj) is given by, 

Q{x- iu) = i x Y,(-l) ax G aa (x; iu), (6.2) 

a 

with a = 0(1) for R(L). As shown in Ref.f/J], 

Q{Q-iuj) = iY,{-mKm\\ (6.3) 

n,m 
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and we shall calculate p(e) by performing analytic continuation iu — > e + iuo. 
The low-order corrections of A in the averaged Green function is given as 

oo oo 

^(0;iw) = <X;(-l) n <o+<E(- 1 ) n <i"- 

n=0 n=0 

n 

+2i £(-l)> n (A0 n , o + 6<f> n , ) + i E(-l)>n,i 

n n 

+0(A 2 ), 



(6.4) 

where each raw corresponds to the order 0(1), O(X), and higher. Inserting the <p n ,m 
and \(fi n ,rn + 5<j>n,m, an d using the formula |TT[ 



/ dn (2(-) $KJ2y/nu) - 2(-)* KMy/nu)) 

JO U> UJ 



1 p/i+l/-l 

4(-)^ +1 ■ — r(i + Ai)r(i + 1/), (6.5) 

ZUJ 1 + fl + V 



we have 



[•OO 

Y^(-l) n 4> 2 n ~ 2cic 2 / dn 0(n, m )(j>(n, -m ) 
„ jo 



u;(l -a; 2m o) 2 



(6.6) 



v~v \»i i / a v x Am 2 cj 2m ° ,2m (l + w 2m °) 2 N . „. 

Vf-iy> 2 ~ 2Am °" 2m ° ( 2m ° 2 ) re a 

V j ^"(1-^0)2^ w 3 J' ^ 



Higher-order terms can be calculated in a similar way. We are interested in the 

averaged Green's function in the limit m — »■ 0, since the derealization transition 

occurs there. Collecting the higher-order terms, 

_ ^ 
0(0; iw) =7Tn \i 
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3' 



+A(-^- 

. 2 , 13i 

+A* 7— ^ 777 + 77"^ ^ + 



-A 



30a;(lna;) 2 3uj(\iiuj) 3 2cj(1iicj) 4 ' 
1 i 



210a;(lnu;) 2 

+0(A 4 ). (6.9) 

Each term in the averaged Green's function has the form ^^V^ . We perform analytic 
continuation in order to obtain the density of states p(e) such as 

Im— -> -Im— — ~ (-1)" : , (6.10) 

w(ma;) a e(lne-if) Q v ; 2e|lne| Q+1 ' v ; 

for small e. Recovering the constant g, the density of states is obtained as follows, 

p( e ) = * 

2t 9 ^t 9 f 

-2g\ j 

24(ln^) 4 

„ 2^2/ 13 1 1 



30^(ln^r 24(ln^) 4 4 (ln£) 



_ 8 ^A 3 7 

210i(lnt) 3 

+0(G?A) 4 ). (6.11) 
The above expression ( |6.11| ) is the main result of this paper. 



7 Discussion 

In the previous sections, we have obtained density of states of fermions interacting 
with quenched disorder of nonlocal correlation. The result is given by (|6.11| ) and Figs.l 



and 2. It is obvious that unit of the energy e is g, the magnitude of disorder fluctuation, 
and p(e) scales with (gX). In the telegraph process of random disorder 0, f and A 
correspond to the mean height and width of fluctuating mass m(x), respectively. 
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Recently we studied the eigenvalue problem 



hip = Eiff (7.1) 

by numerical calculation for telegraphic m(x), and found that all the eigenvalues are 
determined by the combination (gX) in unit ff||12|. 

In Fig.l, we show p(e) for various A with a fixed g. It is easily seen that for 
e smaller than e ~ <?e~«*, p(e) is a increasing function of A whereas for e > e 
it is a decreasing function. That is, states near the band center are enhanced by 
the nonlocal correlation of disorder, which partially suppresses random fluctuation of 
disorder. From the above result, we expect that the number of extended states near 
the band center is increased by the existence of such kind of disorder correlation. 

The above expectation can be verified by solving (|7.1|) . In order to study the field 
equation (|7.1|) in the multi-solitonic background of m(x), transfer-matrix formalism 
is very useful | fL2]| . Using this method, we have obtained solutions to ( [7. 1| ) in various 
configurations of m(x). We think that from these explicit form of solutions we can 
understand characteristic properties of the present random system and their origins, 



e.g., power-law decay of Green's functions, complex multi-fractal scaling|13[, etc. 
Results in [|l^j are summarized as follows; 



1. For quasi-periodic m(x) almost all states are extended and quasi-periodic. Peaks 
of the wave functions appear quasi-periodically. This is a reminiscence of Bloch's 
theorem although we are not imposing the periodic boundary condition on the 
wave function. This result remains correct even for tjiq = (m(x)) ^ 0. 

2. For large random fluctuation of disorder close to the white noise, almost all 
states tend to "localized" even for small or vanishing tjiq. They are divided 
into two categories; the first one simply has one peak or peaks close each other, 
whereas the second one has more than two peaks which are spatially separated 
comparatively large distance. The first one is obviously genuine localized state. 
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3. On suppressing randomness of disorder, the number of states of the second 
category increases. 

We think that the above results support our conclusion from the density of states, i.e., 
the second category is main part of extended states and the number of extended states 
is increased by nonlocal correlation of disorder. We also expect that on increasing 
correlation of disorder, quasi-periodic states in quasi-periodic m(x) come to contribute 
to physical quantities as extended state. 

Sometimes properties of the present disordered system are explained by using the 
single zero-energy wave function ijj± cx e =t ^' ^ x m( - x -* |14]. However as the behavior 



of the density of states shows, a large number of low-energy states are contributing 
to physical quantities. Then it is important to see if these "typical" low-energy 
modes have almost similar behaviors to those of ip±. Details will be reporteded in a 
future publication |I2]. We also calculate Green's functions in the present system [jTU | . 



All these results are quite useful for understanding localization and derealization 
transition. 
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Appendix 



A Higher-order terms 



Higher-order corrections to the averaged Green's function are calculated similarly to 
the correction of 0(X). It is easy to see that the limit, m — > 0, and the integration, 
Eq. fl6.5p , are commutable with each other, and thus we take the limit m — > first. 



From the definition of (j) n ,mi Eq. (|5.15|) , and taking that limit, we have 
* 1 



J n,m 



2an In a; 



(n,m, 0)(1 + 



\n+m\ 



m=l,2,3,... , 



with 



0(n,l,O) = -a V2Xu ■ 2(-) 5 K 1 (2 v ^7), 

UJ 

0(n,2,O) = \/2a \(-uj ■ 2(-)^K 1 (2y/^) + to 2 ■ 2(-)^K 1 (2 v / ^)\ 

V C<J UJ J 

2 3 / 77, l 72 

0(n,3,O) = ^a A2 w • 2(-)2iTi(2v / m^) - 3cj 2 • 2(-)ir 2 (2 v / ^) 

V3 V a; ' uj 



UJ 



Solving Eq.( [5.32|) , we have the higher-order terms in 

1 



A0 n ,m + 

with 

V(n,l,0) = 
V(n,2,0) = 



2a n lna; 



ip(n,m,Q)(l + (-1 



(m = l,2,3, 



V2a Ai f a; • 2(-)* K^y/no) - -w 2 • 2{-)K 2 (2y/ruj 
\ uj 2 uj 



+^.2(-)lK,(2 

6 UJ 



2 

.5 1 i\/nuj 



^a (uj ■ 2(-) 1 *K 1 (2V^) - ^-uj 2 ■ 2{-)K 2 {2^) 
6 \ uj 2 uj 

+5uj 3 • 2(-)5K 3 (2v^J) - -uj 4 ■ 2(-) 2 K 4 (2, 
uj 4 uj 



(A.1) 



(A.2) 



^nu7)j . 
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From the above results, we have the higher corrections to the averaged Green's 
function; Terms in 0(A) are 

2E4,»(A4,o + ^,.)(-D" = + ^3), 

terms in 0(X 2 ) are 

„ \30u;(mu;)V 

£(A</>„, + ^n, ) 2 (-l) n = A 2 (^^7^ + ^^7^ + ^> ~ 7T~ vj) > 

„ V 15 cc;(m c^;)^ 3w(lnw) d 2u(lnw)V 

and terms in (A 3 ) are 

? ^_ ir ._ A 3|L_l_, 

2EU^ + WH)- = A'(^), 

?^-' + *^ 2 (- 1 »"-- A3 (^^)- 

Combining these terms, we have the averaged Green function Eq.( |6.9| ). 
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Figure captions 



Fig. 1 : The density of states, as a function of the energy, whose A is varied from 
0.00 to 0.80 at g — 0.5. The axis of abscissa stands for the energy (e) and the one of 
column for the density of states (p(e)). 

Fig. 2 : The density of states whose g is varied from 0.25 to 2.0 at A = 0.1. 
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